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Abstract. Density fluctuations in the matter distribu- 
tion lead to distortions of the images of distant galaxies 
through weak gravitational lensing effects. This provides 
an efficient probe of the cosmological parameters and of 
the density field. In this article, we investigate the statisti- 
cal properties of the convergence due to weak gravitational 
lensing by non-linear structures (i.e. we consider small an- 
gular windows 9^1'). Previous studies have shown how 
to relate the second and third order moments of the con- 
vergence to those of the density contrast while models 
based on the Press-Schechter prescription provide an esti- 
mate of the tail of P(k). Here we present a method to ob- 
tain an estimate of the full p. d.f. of the convergence k. It is 
based on a realistic description of the density field which 
applies to overdense as well as underdense regions. We 
show that our predictions agree very well with the results 
of N-body simulations for the convergence. This could al- 
low one to derive the cosmological parameters (fi m , J7a) as 
well as the full p. d.f. P(5r) of the density contrast itself in 
the non-linear regime from observations. Hence this gives 
a very powerful tool to constrain scenarios of structure 
formation. 



gravitational lensing 



Key words: cosmology: theory 
large-scale structure of Universe 

1. Introduction 

One important goal of cosmology is to obtain the proper- 
ties of the dark matter density field which eventually led 
to the formation of galaxies, clusters and other astrophys- 
ical objects we observe today. The traditional method to 
obtain some observational constraints on the distribution 
of matter is to build large surveys of galaxies, which map 
the distribution of light, and to use a theoretical model 
to link these galaxies to the density field. However, this 
indirect technique presents the disadvantage to introduce 
significant uncertainties due to the relation between mass 
and light one needs to introduce. Hence it is clearly im- 
portant to build other independent methods which probe 
in a more direct fashion the distribution of matter. Such 



a tool is provided by the gravitational distortion of light 
rays coming from very distant sources. Indeed, the density 
fluctuations along (or close to) the line of sight amplify and 
shear the images of distant galaxies. Then, one can mea- 
sure for instance the shear induced by weak gravitational 
lensing from the observed ellipticities of these galaxies, if 
we assume that the source galaxy ellipticities are uncor- 
rected. Since such effects are purely gravitational they 
probe the total matter content of the universe, dark mat- 
ter as well as baryonic matter. In particular, the statisti- 
cal properties of the convergence or the shear induced by 
gravitational lensing can be directly linked to the charac- 
teristics of the density field. Many authors have already 
considered this problem, focusing especially on the second 
and third order moments in the quasi-linear regime rele- 
vant for large angular windows (e.g., Blandford et al.1991; 
Miralda-Escude 1991; Kaiser 1992; Bernardeau et al.1997; 
van Waerbeke et al.1999). The advantage of such large 
angular windows (9 <; 10') is that one can use rigorous 
perturbative results to describe the properties of the den- 
sity field (e.g., Bernardeau 1994). However, smaller angu- 
lar scales also present a strong interest as they probe the 
non-linear regime which is directly linked to astrophysical 
objects like galaxies, Lyman-a clouds,... Thus, they could 
provide a check on the scenarios used to describe structure 
formation as well as the building of galaxies. Moreover, 
going to smaller scales is also an observational advantage 
since the signal is much easier to measure because it be- 
comes larger while the systematics remain basically the 
same (e.g., van Waerbeke et al.1999). 

Unfortunately, the non-linear regime is much more dif- 
ficult to handle theoretically. Thus, previous theoretical 
studies were restricted to second or third order moments, 
using fits to the results of N-body simulations for the 
non-linear evolution of the power-spectrum (Jain & Seljak 
1997; Hui 1999), or described the density field as a collec- 
tion of virialized halos, obtained from the Press-Schechter 
prescription (Press & Schechter 1974), which allows one to 
model the tail of the distribution of the convergence (e.g., 
Porciani & Madau 1999). However, these works did not 
provide a realistic description of the full probability distri- 
butions of the convergence. Indeed, this requires a consis- 
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tent model for the density field which can describe under- 
dense as well as overdense regions. On the other hand, sev- 
eral numerical studies have been performed using N-body 
simulations (e.g., Jain et al.1999). In this article, extend- 
ing an earlier work relevant for point sources (Valageas 
1999a), we present an analytic method to obtain the p.d.f. 
of the convergence. This is based on a scaling model devel- 
oped in Balian & Schaeffer (1989) (see also Bernardeau & 
Schaeffer 1992; Valageas & Schaeffer 1997) which has been 
seen to provide a good description of the non-linear den- 
sity field through comparisons with numerical simulations 
(e.g., Valageas et al.2000; Bouchet et al.1991; Munshi et 
al.1999). 

This article is organized as follows. In Sect]! we briefly 
present the model we use for the density field and we in- 
troduce a few statistical tools. Next, in Sect.|| we describe 
how we obtain the p.d.f. of the convergence k, smoothed on 
small angular scales (9 1'). We consider the case where 
all sources are at the same redshift Zg clS well as a broad 
redshift distribution. Then, in Sect. |] we describe our nu- 
merical predictions, which we compare to available results 
from N-body simulations, for three cosmologies. Finally, 
in Sect.|5| we briefly consider the dependence of the weak 
gravitational lensing effects on the redshift of the sources. 

2. Properties of the density field 

In order to obtain the characteristics of the weak gravita- 
tional lensing effects we need a detailed description of the 
density field which gives rise to these distortions. Hence 
we briefly recall here the formalism we use to describe 
the density fluctuations. See Balian & Schaeffer (1989) for 
the details of our approach (see also Valageas & Schaeffer 
1997; Valageas et al.2000). We shall use the same methods 
to obtain the properties of the weak gravitational lensing 
effects. 

Our main tool to obtain the probability distributions 
of the various quantities we consider will be the generat- 
ing functions defined from the moments or the cumulants 
of these variables, which are widely used in statistics (see 
also Balian & Schaeffer 1989). Thus, the probability dis- 
tribution P(fJ,) of a random variable fi can be derived from 
the generating function ip^ (y) (similar to a Fourier trans- 
form) defined from the moments (fi p ) (provided they are 
finite) by: 



Mv) = J2 

p=0 



(-l) p 



since we have: 

r +ioc dv 

J-ioo 2m 



(1) 



(2) 



Moreover, it is often useful to introduce the generating 
function ^(y) defined from the cumulants (fJ- p ) c by: 



<*v(y) 



p 



(3) 



which verifies: 



%(y) = hx[i/> lt {y)] 



(4) 



In the following, we consider variables such that their 
mean is zero. Then it is convenient to define a new gener- 
ating function (p^{y) by: 

- ( ^ f HTjM!k v i x n M 

Mv) = l^— ] — (At 2) P -i y > o*) = ° ( 5 ) 

Note that in this case (/i 2 ) = (p, 2 ) c since (/i) = 0. Then, 
using (0) and @ we obtain: 



$ (y) = e -Mv£»)/t» with ^ = (A* 2 ) 



and: 



-Hoc 



dy 



,[»v-0r(y)]/£» 



(6) 



(7) 



This formalism allows us to determine for instance the 
probability distribution P{b~n) of the density contrast 8r 
within spherical cells of volume V, radius R: 



= / — s(r) with 5 = P -^- 
Iv V p 



(8) 



in terms of the many-body correlation functions 
£p(ri) T p ) (here ~p is the mean density of the universe). 
Indeed, from the definition of the correlation functions 
(e.g., Peebles 1980) we can write the cumulants (S^) c as: 



P>2 : (a# c = £p 



d 3 ri...d 3 r p 
VP 



€ p (ti,...,t p ) (9) 



Next, we define the parameters S p and the generating 
function <p(y) which corresponds to (^) by: 

= E -^r~ s p y p > ^ = =^1.^ = 1 (io) 



p=2 



£2 



Then, from (^) one obtains (White 1979; Balian & Scha- 
effer 1989): 



-Hog 



P( Sr ) = I -% e \s*v-mvt 

2nit; 



(11) 



where we note the mean correlation £ 2 at scale R as £. 
The interest of the generating function <p(y) is that in the 
quasi-linear regime the parameters S p remain finite in the 
limit £ — > so that one can estimate P(6r) from (11) 



using the function (p(y) obtained by rigorous perturba- 
tive methods for £ — > (Bernardeau 1994). Moreover, in 
the non-linear regime which we consider in this article the 
stable-clustering ansatz (Peebles 1980) also implies that 
the parameters S p reach a finite limit in the highly non- 
linear regime £ ^> 1 for an initial power-spectrum which 
is a power-law (e.g., Valageas & Schaeffer 1997). In this 
case, once ip(y) is measured at one scale and time in the 
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highly non-linear regime, the density probability distribu- 
tion P(Sr) can be obtained for any time and scale in the 
non-linear regime provided that one knows the behaviour 
of £. Note that the stable-clustering ansatz provides the 
behaviour of the correlation functions in the highly non- 
linear regime since it means that at these scales one has 
(Peebles 1980): 

^(Xr 1 ,...,Xr p ;a) = a 3 ^ A^" 1 ) | p (n, r p ) (12) 

where a(t) is the scale- factor and 7 is the (local) slope of 
the two-point correlation function. Note that for an initial 
linear power-spectrum which is a power-law P(k) oc k n we 
have if stable-clustering is valid: 



7 



3(3 + n) 



(13) 



The consequences of the scaling-laws (|12|) were studied in 
details in Balian & Schaeffer (1989). In practice, one does 
not directly measure <p(y) itself from numerical simula- 
tions or observations but the p.d.f. P(6r) from which it 
can be derived using (f[l|). More precisely, in the highly 
non-linear regime one considers the variable x and the 
generating function tp(y) defined by: 



•'• = 1 + - Sr , <p(y) = v + <p{y) 



(14) 



Then, using ([H]), for sufficiently "large" density contrasts 
the p.d.f. P{Sr) can be written as (Balian & Schaeffer 
1989): 

I » 1 , (1 + S R ) » ? : P(6 R ) = Jj h(x) (15) 

where the scaling function h(x) is the inverse Laplace 
transform of <p(y): 



h(x) = - 



which obeys: 



P>1 ■■ s p 



x p h(x) dx 



S\ — S2 



(16) 



(17) 



In (|T^) the exponent uj comes from the behaviour of tp(y) 
at large y, see (A..1). Note that for large £ the range where 
( |l5| ) is valid extends from very low densities (and negative 
Sr) up to infinity. Thus, one can obtain h(x) from the 
counts-in-cells statistics measured in simulations. Then, 
one can derive f(y) from h(x) since (Era) can be inverted 



as: 



<p(y) = 



h(x) dx 



(18) 



This also gives <p(y) from (|l4|). We recall in App.A the 
behaviour of P(8r) obtained for simple models for h{x) 
and ip(y) which are consistent with numerical results. The 



function h{x) has been measured in the non-linear regime 
for various power-spectra (CDM models and power-laws) 
by several authors (Valageas et al.2000; Bouchet et 
al.1991; Colombi et al.1997; Munshi et al.1999). In partic- 
ular, although Colombi et al.(1996) found a small scale- 
dependence other authors found that the numerical results 
were consistent with h{x) being scale-independent in the 
non-linear regime. Thus, in the following we shall use the 
scaling function h{x) obtained by Valageas et al.(2000) for 
n = —2. Note that h{x) depends on the power-spectrum 
and it must be obtained from numerical simulations since 
there is no known method to derive analytically h(x) (see 
App.A). In our case, we choose the scaling function mea- 
sured for n = —2 because as we shall see in Sect. 4.2 and in 
Fig.H most of the contributions to the weak lensing effects 
we study in this article come from scales where n ~ — 2. 
As described in Valageas (1999b) this model provides the 
simplest realistic description of the non-linear density field 
which takes into account substructures and underdense as 
well as overdense regions. 

Then, to obtain the properties of the non-linear den- 
sity field we only need to model the evolution of the 
two-point correlation function z) (or the non-linear 
power-spectrum) from the initial linear power-spectrum. 
To this order we use the fits given by Peacock & Dodds 
(1996) which give the non-linear power-spectrum P(k) 
from its linear counterpart. Note that the evolution of 
£ measured in numerical simulations (e.g., Valageas et 
al.2000; Peacock & Dodds 1996) is consistent with the 
stable-clustering ansatz. 

3. Distortions induced by weak gravitational 
lensing 

3.1. Shear tensor 

The gravitational lensing effects produced by density fluc- 
tuations along the trajectory of a photon lead to an appar- 
ent displacement of the source and to a distortion of the 
image. Thus, light coming from a direction 9 is deflected 
by a small angle 59. However, the observable quantities 
are not the displacements 59 themselves but the distor- 
tions induced by these deflections, which are given by the 
symmetric shear matrix (e.g., Jain & Seljak 1997): 



d5dj 
~d0~ 



= -2 



dx ¥T( — ^ VjVj0(x) (19) 



Here x is the radial comoving coordinate (and \s corre- 
sponds to the redshift z s of the source): 



dx = 



-TT- dz 



^/n A + (1 - n m - n A )(i + z y + a„(i + *) 3 

while the angular distance T> is defined by: 
c/H 



(20) 



V(z) 



sinh (yi - n m - r» A (21) 
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The gravitational potential <f> is related to the density fluc- 
tuations 8 by Poisson's equation: 



(22) 



In ( |l9| ) we used the weak lensing approximation: the 
derivatives ViVj0(x) of the gravitational potential are 
computed along the unperturbed trajectory of the photon. 
This assumes that the components of the shear tensor are 
small but the density fluctuations 5 can be large (Kaiser 
1992). The shear tensor <E>jj is usually decomposed into 
its trace k and the shear components 71, 72, defined by: 



$1,1 + $ 



2.2 



and 



71 



$1 



72 



(23) 



-$1,2 , 7 = 7i + i 72 (24) 



Moreover, the magnification fj, of the source is given by: 
1 



(25) 



(l-«)a_| 7 |a 

so that for small values of the shear components we have: 



K < 1 



1 + 2k 



(26) 



Thus the goal of this article is to obtain the probability 
distribution of the component k of the shear tensor, as de- 
fined in ( [l9| ) . This applies to point sources like supernovae 
or weak gravitational lensing effects which are filtered on 
small angular scales (9 & 1') since we only consider here 
the non-linear regime. 

3.2. Convergence n 

The probability distribution P(p) of the magnification /j, 
as defined in (|2q ) without smoothing was already obtained 
in Valageas (1999a). This directly gives the probability 
distribution P(n) of the convergence k. Here we extend 
this method to the case of a finite smoothing angle 9. This 
allows us to check the validity of our treatment through 
comparisons with available results from N-body simula- 
tions. 

Using (|l^) and ( |22] ) one can show (Bernardeau et 
al.1997; Kaiser 1998) that the convergence along a given 
line of sight is: 



with: 



dx w(x,Xs) Hx) 



Hi V(x)V( X s - X) 



V(Xs) 



(1+*) 



(27) 



(28) 



where z corresponds to the radial distance x- Thus the 
convergence k can be expressed in a very simple fashion as 
a function of the density field. First, we can see from ( p7| ) 



that there is a minimum value K m in{ z s) for the conver- 
gence of a source located at redshift z s , which corresponds 
to an "empty" beam between the source and the observer 
(5 = — 1 everywhere along the line of sight): 

« mi „ = -^F s ( Xs ) (29) 



with 

F s (Xs) 



dxw(x,Xs) (30) 







Next, we define the "normalized" convergence k by: 



(31) 



which obeys k > —1. If one smoothes the observations 
with a top-hat window in real space of small angular radius 
9 one rather considers the mean "normalized" convergence 
kg: 

Here C is a vector in the plane perpendicular to the line 
of sight (we restrict ourselves to small angular windows) 
over which we integrate within the disk \C,\ < 9 (we note 
this by the short notation J Q ) . Thus x is the radial coordi- 
nate while T> C is the two-dimensional vector of transverse 
coordinates. 

3.3. Probability distribution P(k) 

Next, in order to obtain the p.d.f. P(kg) of the conver- 
gence we simply need to derive the cumulants (kg) c . This 
will provide the parameters Sk b , p and the generating func- 
tion ipgAy) in a fashion similar to ( ^p| ) and (jl^). Indeed, 
from (32) we have (kg) = 0. As noticed by several au- 
thors (e.g., Bernardeau et al.1997; Hui 1999; Munshi & 
Coles 1999) the cumulants ( be simply expressed 

in terms of the p— point correlation functions of the den- 
sity field through (p2|). Thus, we write: 



e p 



J- J- 7T0 2 



n 



1 



(33) 



w{Xi,Xs 



8(Xi>Vi Ci) 



From the definition of the correlation functions (Peebles 
1980): 



£ p (xi, ...,x p ) = (6 (xi) ...S (xp)> c 
we obtain: 



(34) 



Xa , w(xi,Xs) f Xs Xi tj, w(xi+Xi,Xs) 
d X i p / [[dXi p 

r s J-Xl i=2 s 



' P d 2 d 



Xi 



Xi +x P 
v P C P 
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where we made the change of variables Xi ~ * Xi + Xi 
for i > 2. Since the correlation length (beyond which the 
many-body correlation functions are negligible) is much 
smaller than the Hubble scale c/H(z) (where H(z) is the 
Hubble constant at redshift z) only values of Xi which 
obey \xi\ -C c/H(z) contribute to the integral over x% (for 
i > 2). Hence we have w(xi+Xi,Xs) - w(xi>Xs), A - ^i 
and we can push the integration boundaries over Xi (Jot 
i > 2) to infinity. Thus, we get: 



w(x,Xs) 
F s 



P r oo P 







(=1 



1 i=2 



Xp 

T>Ci""' 



(35) 



where we used the fact that £ p is invariant through the 
translation (— XijO) over the p points x^: x% ~* Xi ~ Xi- 
Although the points (xi,P O cover a cylinder rather than 
a sphere, we approximate the integral over the p-point 
correlation function by: 



r p r 9 p d 2 t- 



s P IT 1 



(36) 



in a fashion similar to £ p = S p £; 
where we defined: 



-p-i 



, see © and ©, 



t:9 2 t:9 2 ^\VCi'VC 2 ,Z 



(37) 



We discuss in more details in App.B the approximation 
( |36| ) for the case of a tree-model for the p— point correla- 
tion functions. However, our approach only assumes that 
the stable-clustering ansatz is valid (i.e. the coefficients S p 
are constant with time at a given physical scale). In par- 
ticular, the approximation ( |36|) should provide reasonable 
results even if the p— point correlation functions are not 
exactly given by a tree-model. Thus, using d36| ) we write 
the cumulants (B5|) as: 



(38) 



dX ( f ) S p P- 1 



From the definition (||) and (|3^) we obtain the generating 
function (pk e {u)- 



<p* e {y) 



dx — v [yyj- 



where we introduced the variance: 



(39) 



(40) 



Throughout this article, ip (as in the r.h.s. of ([39|)) and 
tp, without subscript, refer to the generating functions de- 
fined for the density field in dj) and App.A. 



Note that the variance ^ g we obtain in ( p0| ) docs not 
rely on the approximation ([36]) and it is exact (within the 
weak lensing approximation). We only use (|3^) to get an 
approximation for the higher-order parameters Sf- 9jP with 
p > 3. From (]?]) we obtain the p.d.f. P(kg) and P(ng) as: 



P(k e ) = 



and 



P(k 8 ) 



-\-ioo 



dy 



P(He) 



,[KeV-<Pk e (v)]/in B 



(41) 



(42) 



The relations ([id]) and ([39]) are the main results of this 
article. They allow us to derive the properties of the con- 
vergence from the counts -in- cells statistics in a straightfor- 
ward fashion. The term I K defined in ( ]37| ) can be written 
in a more convenient form as a function of the power- 
spectrum P(k) defined by: 



<<5(k 1 )5(k 2 )) = P(fci)5 D (k 1 + k 2 ) 
Thus, using: 



(43) 



£ 2 (x) = J d 3 ke lk x P(k) , k.x = fc|| X + k ± .I?C (44) 

where fcii is the component of k parallel to the line of 
sight while kj^ is the two-dimensional vector formed by 
the components of k perpendicular to the line of sight, we 
obtain 



dk^l w2{me) 
/o k k 

where we defined at redshift z: 
A 2 (k;z) = 4vrfc 3 P(fc;z) 
and: 

2 



W (Vk9) 



Vk6 



Ji (Vk9) 



(45) 



(46) 



(47) 



Here J\ is the Bessel function of the first kind of order 1 . 
Note that throughout this article x and k refer to comov- 
ing quantities. 

From ([38]) we can also obtain the coefficients S Ke lP de- 
fined in a fashion similar to (]To|): 



*0/C 



i2— p v^e/c 



(48) 



In particular, we get for the skewness S Kg ^ the expression: 



s K = — s 3 f dx w3l \ 

[J dx W 2 I K ] ' 



(49) 



We can see from ( ]l9| ) that S Kgl 3 is (almost) independent 
of the normalization of the power-spectrum, while it shows 
a strong dependence on the cosmological parameter f2 m . 
Thus it could be used to measure Q m , see Bernardeau et 
al.(1997) for such a study in the quasi-linear regime. 
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It is possible to get some important information about 
the generating function (p kg (y) by direct inspection of (|39|). 
First, we note that if we define the function ip Kg (y) by: 



V£ e (y) = y + <Pk e (y) 



(50) 



as in (jlil), then from ( |39| ) we see that <Ph g (y) and ip(y) 
bear the same relation as <^« e (y) and (p(y): 



ax — p [y- 



(51) 



Moreover, using (gl|) the p.d.f. P(kg) can be expressed as 



P{kg) 



+200 



dy 



,[(l+ft8)»-(P* e (i/)]/£* e 



(52) 



In this article, we use ( |39| ) and (|4l|), rather than (51) and 
(|52]) , because for small values of the variance ^„ the p. d.f. 
P(kg) still looks like a gaussian (except in the tails of 
the distribution). This is most clearly seen in (41) since 
'fUgiy) — — y 2 /2 + .... Moreover, the use of (p(y) ensures 
that there is no linear term in y in the exponent in (]4l]) . In- 
deed, a small numerical error in the linear term of ip Kg (y) 
(which must satisfy p Kg {y) = y — y 2 /2 + ■•■) would show 
up as a non-zero mean (kg) if one uses (|52|). However, for 
theoretical considerations ( pi"]) a nd ( [52] ) are more conve- 
nient. Thus, from ( |5l| ) and (A.l) we see that <p Kg (y) and 
<p(y) have the same power-law behaviour for large y: 



R-e(y) -> +00 : ip kg (y) ~ a kg y 1 



a K „ > 



(53) 



In our case uj ~ 0.3. Here we also used the fact that I K 
is always positive, as shown by (^BJ). Then, from (53) and 
(|52]) we check that P(kg) = for Kg < —1. Indeed, for 
Ke < -1 we can push the integration path in (|52| ) towards 
the right, Re(j/) — ■> +00, so that the integral vanishes. 
Thus, our approximation ^b\ ) has preserved the fact that 
P(kb) = for Kg < K min . Moreover, using ( j53"l ) in (|52] ) we 
obtain the behaviour of the cutoff of P(n§) for Kg — ► K m i n : 



P(«e) ~ exp 



w a 



l/a 



k 9 



(1 - w)|K m in| 



-(l-w)/a 



(54) 



where we did not write multiplicative power-law factors. 
Next, using (A. 2) we see from ( ]39j ) or ( |5l| ) that (pa g {y) 
and (pa e (y) have a singularity (branch cut) at y s Kg given 
by: 



lis 



Us 



max 

x 



< 



(55) 



The location of this singularity is important since in the 
numerical integration of ( fH] ) one must make sure that the 
integration path does not cross the branch cut (i.e. the 
integration path crosses the real axis at a point y such 
that y > y s ,k g )- Moreover, as can be seen from (|d|) the 



existence of the singularity at y s ,k g implies that P(kg) 
shows an exponential tail for large kg : 



kg > 1 : P(kg) ~ e -^».* 9 l 

so that P(Kg) obeys: 



with: 

ng.s = 



£k 8 



3fL 



\ys,a e \ 



x s max(w/ K ) 

x 



(56) 
(57) 

(58) 



Here we used x s = l/\y s \ as in (A. 3). Note that the ex- 
ponential cutoff of the p.d.f. P(ng) is stronger for low- 
density universes and a lower normalization of the power- 
spectrum P(k) (through I K ) since Kg s gets smaller. Of 
course, this is due to the fact that in such cases there are 
fewer high-density massive objects (which are the ones 



which can produce a high convergence Kg). Note that (58) 



shows that the cutoff of the p.d.f. P(Kg) strongly depends 
on the cutoff a; s of the density contrast distribution, so 
that a measure of the tail of P(k$) provides a direct esti- 
mate of the tail of P{Sr) (but see the discussion in Sect.||). 

The case without smoothing is recovered in all previ- 
ous expressions by taking the limit 9 — > 0. In particular, 
( |39|) and (^lj) are unchanged while in (^5|) the factor W 2 
is equal to 1. Of course, we could also perform the same 
calculation directly without smoothing, so that k only in- 
volves the integration along one line of sight, as in (p7|), 
and we would recover these results. 

3.4- A convenient approximation 

We can check from (|39| ) that the expansion about y = of 
the generating function (p Ke (y) satisfies: (p Kg (y) = —y 2 /2 + 
as implied by the definition (||). Moreover, from (^o|), 
( |40l) and ( [39|) we see that (pa g (y) should be close to <p(y) 
since the factor (wI K )/(F s £ Kg ) has typical values of order 
unity. Thus, a simple approximation is to use: 



<pa e (y) - <p(y) 



(59) 



which simplifies somewhat the numerical calculations. 
Note that the approximation ( |59| ) actually means that we 
use S Kg ^p — 5*p, i.e.: 

^K g ,p — | K min| P Sp (60) 

The approximation ( |59| ) also implies that P(ng) shows an 
exponential cutoff for large positive kg with a parameter 
similar to Kg tS given by: 

3fL 



Kg.. 



•^s (Ps£,fc g ) 



(61) 



In fact, using (|1 1|) , the approximation ([59|) means that 
the p.d.f. P(kg) is directly given by the p.d.f. P(5r) of 
the density contrast as: 



P(ko) = 



1 



P 6 



Kg 



(62) 



P. Valageas: Statistical properties of the convergence 



7 



This clearly shows that the p.d.f. P(kg) is an efficient tool 
to obtain the properties of the density field. Indeed, we 
shall check in Sect. 4.4 that the approximation ( p9| ) gives 
very good results and the relation (|6^) provides a straight- 



forward method to derive P{5r) from P(ke). Note how- 
ever that even for ^ e < 1 the p.d.f. P(5r) used in ( |62[ ) 
corresponds to the non- linear regime (for 9 ^ 1'). Hence 
it is not the p. d.f. of the density contrast measured at the 
time when £ = since the values of the coefficients S p 
in the linear and non-linear regimes are different. 



We shall see in Sect. 4.4 that the approximation (p9[) 



gives very good results. In particular, we noticed above 
in ( J53| ) that (p^ g (y) and (p{y) obey the same power-law 
behaviour for large y. However, the location of their sin- 
gularity along the real axis is slightly different, see (|55|), 
and the exponent uj s of the power-law prefactor at large x 
of their associated inverse Laplace transform h(x), defined 
as in (|i"6|), differs by a factor 1/2, as discussed in Valageas 
(1999a) (see the system (43)). More precisely, we obtain 
from fl55|): 



\Vs\ _ J d\w max(w/ K ) 



I d X 



I K 



> 1 



(63) 



which means that the singularity y s ,k e is closer to than 
y s . This implies that the exponential cutoff of P{kq) is 
slightly smoother for the more accurate expression ( p9[ ) 
than for the approximation (|59|), see (58). Next, we noticed 
above that ( |59| ) implies that the skewness Sk e ,3 of the 
normalized convergence kg is given by Sa Bt 3 = S3 while 
the expression (E9) gives: 



s 3 



Jd\w Jdxw 3 I 2 
[J dx w 2 I K ] 



> S 3 



(64) 



for the more accurate calculation (39). The inequality is 
obtained from the Cauchy-Schwarz inequality applied to 
the scalar product {f\g} denned by: (f\g) = J dx f.g for 
the functions / = w 1 ^ 2 and g — w 3 ^ 2 I K . Hence the skew- 



ness of the convergence given by the approximation ( |59| ) is 
slightly too small. The inequalities (|6^ ) and ([34|) translate 
the fact that the approximation ( |59| ) corresponds to an 
average along the line of sight of the kernel {wI K )/{F s ^ g ) 
which appears in (^), hence to an average over the var- 
ious p.d.f. of the weak lensing effect arising from differ- 
ent planes along the line of sight. On the other hand, 
the expression ([59]) takes into account the variations of 
this kernel so that the p.d.f. of the total convergence is 
given by the exact convolution of the various p.d.f. associ- 
ated to the elementary contributions due to the successive 
mass planes along the line of sight (which are discretized 
in numerical ray-tracing simulations). Then, the tails of 
the p.d.f. P(ng) are sensitive to the highest tails among 
these successive p.d.f, as shown by the maximum which 
appears in (|55|). This implies that the tails of P(ng) given 
by ( |39| ) are smoother than the cutoffs implied by the ap- 
proximation (|59]). This leads to a smaller cutoff Kg s and a 



smaller skewness for this approximation, as shown in (63) 
and @. 

3. 5. Redshift distribution of the sources 

In the previous sections we considered the case where all 
the sources are located at the same redshift z s . This is 
convenient to compare our predictions with the results 
of numerical simulations and to study the dependence on 
redshift of the weak gravitational lensing effects. However, 
in practice one observes galaxies over a finite range of red- 
shifts in order to get a sufficiently large number of sources 
to perform a statistical analysis. Hence we show in this 
section how our results can be applied to a "broad" red- 
shift distribution of sources n(z s )dz s . Here we normalized 
n(z s ) to unity: 



dz s n(z s 



1 



(65) 



where z max is the maximum redshift of the sources. Then 
the mean convergence 77 observed in one direction on the 
sky is the average of the convergence k defined in ( |27| ) over 
the redshift distribution of the sources: 



dz s n(z s ) k(z s ) 



(66) 



which gives: 



3ft n 



dz s n(z s ) 



Xs{Zs) 



dx w{x,Xs) 5(x) (67) 



The expression (BTl) can also be written: 



30, 



with: 

w(x) ■ 



*(x) 



dx w(x) 8(x) 



dz s n(z s ) w (x,Xs(z s )) 



(68) 



(69) 



Thus, we obtain exactly the same expression as in (27) 
with Xs replaced by x max and w(x, Xs) by w(x). Hence all 
our previous results remain valid, after we perform these 
two changes and F s is replaced by F s (compare with ( |30| 
obtained from: 



F H = 



dx w(x) 



(70) 



In particular, the relations ( pl| ) and ( |39| ) are still correct 
after we make these straightforward changes. Hence our 
results also apply to the p.d.f. P(~R~e) of the convergence 
observed from a broad redshift distribution of sources. 



Moreover, we noticed in Sect. 3.4 that a convenient approx- 



imation to the generating function if^ e (y) (defined for a 
redshift z s of the sources) is simply ip(y), see (|59|). This 
is checked below in Fig.|[ In the case of a broad redshift 
distribution of the sources we can still perform the same 
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approximation (|39|) since we noticed that the relation ( p9[ ) 
still holds (with the changes of notations described above). 
In fact, this could be expected since <p(y) is independent 
of z s . This means that, to a good approximation, the full 
p.d.f. P(jZq) of the convergence exhibits a specific scaling 
property. Indeed, it is described by a unique generating 
function (p(y) (which characterizes the properties of the 
underlying density field) which is independent of the cos- 
mological parameters, of the angular scale 9 and of the 
redshift distribution of the sources, and by two numbers 
7t m i n and £k 8 (which describe the lower bound and the 
variance of this p.d.f.) which contain all the dependence 
on the cosmology, on 9 and on the redshift distribution 
of the sources. This holds as long as one probes the non- 
linear regime up to z max (i.e. 9 must be small enough). In 
particular, the two p.d.f. P(kg) of the normalized conver- 
gence kg obtained for two different rcdshifts z s \, z S 2, of 
the sources and two angular windows 0\, 92, chosen such 
that the normalized variances are the same, = £/j 92 , 
should nearly superpose. 

Note that in the treatment presented in this section 
the redshift distribution of the sources is described by a 
uniform background population n(z s ). In fact the quantity 
n(z s ) will fluctuate from one line of sight to another and 
these variations Sn should be correlated with the density 
fluctuations along the line of sight, because the galaxies 
are expected to be correlated with the dark matter density 
field. As a consequence, in the expression ( |S8| ) the prod- 
uct n(z s ).S(x) contains a term Sn(z s ).S(x) which should 
be taken into account. However, since the kernel w(x,Xs) 
vanishes for x. = Xs most of the weak lensing effects which 
distort the image of a source arise from density fluctua- 
tions located at cosmological distances (~ c/H(z)) along 
the line of sight from this source. This supresses somewhat 
the contribution of the coupling between the fluctuations 
of the sources 5n(z s ) and of the density field S(x). Then, 
the contributions due to the fluctuations of the sources 
only, which lead to terms like (Sn 2 ), become relatively 
small if the number of sources is sufficiently large. How- 
ever, these effects would certainly deserve a more detailed 
study but this is beyond the scope of this article. 

4. Numerical results 

We now describe the numerical results we get from the 
analytical tools detailed in the previous sections. In order 
to compare our predictions with the available results ob- 
tained from ray-tracing through N-body simulations we 
consider the three cosmological models defined in Tab.|l|: 
a standard CDM scenario (SCDM), a low-density open 
universe (OCDM) and a low-density flat universe with a 
non-zero cosmological constant (ACDM). Here T is the 
usual shape parameter of the power-spectrum. We use the 
fit given by Bardeen et al.(1986) for P(k). We consider 
the weak lensing distortions which affect a population of 
sources at redshift z s — 1. 



Table 1. Cosmological models 





SCDM 


OCDM 


ACDM 




1 


0.3 


0.3 










0.7 


Ho [km/s/Mpc] 


50 


70 


70 




0.6 


0.85 


0.9 


r 


0.5 


0.21 


0.21 



4-.1. Magnitude of the convergence 

First, we consider the rms convergence v/£« e = 
I K min | v^fte" we obtain for various angular windows. Our 
results are displayed in Fig||. 

We see that our results agree rather well with the N- 
body simulations (the points correspond to the average 
over several realizations while the error bars give the over- 
all scatter, see Munshi & Jain 2000). This is expected 
since our prediction for the variance t; Kg only relies on 
the weak lensing approximation (Kaiser 1992) and on the 
prescription given by Peacock & Dodds (1996) to obtain 
the non-linear evolution of the power-spectrum. Since the 
latter is a fit to other N-body simulations, the agreement 
only shows that both sets of simulations are consistent. 
The deviations at small angular scales are due to the fi- 
nite resolution, see Jain et al.(1999) for details. Although 
there remains a small discrepancy for the ACDM scenario 
it is certainly of the same order as the accuracy of the sim- 
ulations. The main feature of Fig.|l] is that the variance is 
smaller for low-density universes, despite their larger nor- 
malization erg of the power-spectrum. Indeed, from (29) 
and ( |40| ) we can write the variance as: 



V 2 



dx w 2 n 



dk A 2 (k) 



w 2 (me) (7i) 



and the figure shows that the strong dependence on f2 m 
is more important than the variation with a$ (for the cos- 
mologies we consider here, which have reasonable param- 
eters). Of course, the variance is smaller for low-density 
universes (at fixed A(fc)) because in such cases there is 
less matter to produce weak lensing distortions. Since the 
difference between an "empty" beam (S = — 1 everywhere 
along the line of sight) and the mean (5 = 0) is also smaller 
- it is proportional to f2 ro as seen in ( |29| ) - the upper 
bound |ft m in| of the magnitude of negative deviations is 
lower than for the critical density universe. We can see 
from Fig|l] that an angular window 9=1' still removes 
some small scale power as compared to the case without 
smoothing (horizontal dashed lines). In particular, the lat- 
ter, which corresponds to the observation of point sources 
like SNela (see Valageas 1999a), is still somewhat beyond 
the resolution of N-body simulations. 
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Fig. 1. The magnitude of the convergence Ke for a source 
at redshift z s = 1 for various angular windows 9 and three 
cosmologies. The solid lines show the rms convergence 
\/^K, g , the horizontal dashed lines show the rms conver- 
gence y/^no obtained without smoothing (9 = 0) and the 
horizontal dotted line in the middle panel shows |K m i n |- 
For both flat universes |K m i n | is larger than 0.055 (0.12 for 
the critical density universe and 0.064 for the model with 
cosmological constant). The data points are the results of 
numerical simulations from Munshi & Jain (2000). 

4-2. Dependence on scale of the contribution to weak 
gravitational lensing 

Next, in order to distinguish the contributions to weak 
gravitational lensing effects arising from different comov- 
ing scales we define the quantity: 



x ° , 2 f°° dk A 2 (k) 

dxW 7T / — — 

o Jo K K 



(72) 



which measures the contribution from the comoving 
wavenumber k to the variance of the convergence £ Ke , see 
( pp| ) and p^). Thus, Aj^(fc) integrated over dlnfc is nor- 
malized to unity for the case without smoothing (9 = 0). 
Hence A 2 e (k) also shows the influence of the finite window 
9. 

We present in Fig.^ our results for three different an- 
gular windows, 9 = 0.1', 9 = 1' and 9 = 10' (shown by 
the dashed curves), as well as the case without smooth- 
ing (upper solid curve). As we can clearly see in the fig- 
ure, the effect of the finite window is to cut the contribu- 
tion from small scales (smaller than the size of the win- 
dow), i.e. from high comoving wavenumbers k. Thus, even 
the smallest window 9 = 0.1' cuts a significant part of 
the small-scale power which contributes to the case with- 
out smoothing. The latter is relevant for point sources, 
e.g. supernovae (see Valageas 1999a for a detailed study 
and the implications for the measure of the cosmologi- 
cal parameters). Since present numerical simulations al- 
ready miss some power for 9 = 0.1' (Jain et al.1999) 
they would significantly underestimate the effect associ- 
ated with point sources (in particular note the extended 
tail at high k). The lower solid lines in the figure corre- 
spond to the case without smoothing where we use the 
linear power-spectrum Ph(k) instead of the actual power- 
spectrum P(k). Thus, we sec from the figure that for small 
angular windows 6 & V the non-linear evolution of the 
density fluctuations plays an important role and must be 
taken into account, as was already noticed by Jain & Sel- 
jak (1997). This shows that the formalism we developed 
in Sect|| applies to 9 ^ 1' where the relevant scales are 
within the non-linear regime. For very large scales 9 <; 10' 
we could also use the same formalism but the parameters 
S p would be different (see Colombi et al.1997). The in- 
termediate regime is rather difficult to model accurately 
since it corresponds to a transitory regime where the co- 
efficients S p evolve (but see Hui 1999). 

For 9 = the characteristic scale which dominates the 
weak lensing distortions (shown by the maximum of the 
upper solid curve in the figure) is set by the shape of the 
initial power-spectrum itself. Thus, from (|72] ) we see at 
once that this corresponds to the wavenumber k m where 
the local slope of A 2 (fc) is 1, i.e. A 2 (fc) cx k. From the 
definition (^) this means P(k) oc k~ 2 , i.e. n = —2 or 
£(r) oc r _1 . In fact, we can see from ( |l3| ) that an ini- 
tial power-spectrum Pi(fc) oc fc~ 2 keeps the same slope 
n = —2 in the highly non-linear regime. However, we see 
from the figure that A 2 o (k) peaks at a value k m somewhat 
larger than the wavenumber fc_2 where the local slope of 
the initial linear power-spectrum Pt(fc) is equal to —2, 
which corresponds to the peak of the lower solid curve. 
This is due to the form of the detailed evolution from the 
linear regime to the highly non-linear regime. However, we 
can check in the figure that for 9 10' the main contribu- 
tion to the weak lensing distortions comes from wavenum- 
bers which are close to fc_2- This justifies the use of ( A. 4 ) 
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Fig. 2. The contribution of various comoving wavenum- 
bers k to the convergence Kg for a source at redshift z s = 1 
for three cosmologies. The upper solid lines show A^ fl (k) 
for the case without smoothing (6 = 0). The dashed lines 
show the cases 9 = 0.1', = 1' and 9 = 10'. A larger angu- 
lar window 9 corresponds to a curve with a smaller upper 
cutoff in k and a lower maximum. The lower solid curves 
show A^ o (k) (i.e. no smoothing) when the power-spectrum 
P(k) is always taken to be the linear extrapolation. 

which gives the parameters S p for the case P(k) oc k~ 2 . 
In fact, the curvature of the actual power-spectrum might 
lead to a small modification of these coefficients S p but 
our approximation has the advantage of the simplicity 
and should provide a reasonable prescription. In order to 
improve our treatment, one would need to measure the 
function ip(y) from the statistics of the counts-in-cells, as 



explained in Sect]^, on the relevant scale for each cosmo- 
logical scenario. 

4- 3. Skewness 

A convenient measure of the non-gaussianity of a peculiar 
p.d.f. is the skewness S3. In particular, we can consider 
the skewness S Ke 3 of the convergence which we obtained 



in Sect. 3.2 



We compare in Fig.|| our predictions to the results 
of numerical simulations from Jain et al.(1999) for three 
cosmologies. Although the numerical results probe scales 
which are not entirely within the non-linear regime (9 <; 
1'), as shown by the variation with 9 of S K0 3 which does 
not appear to have reached its asymptotic value at 9 = 1' 
yet, we can check that they roughly agree with our predic- 
tions. In particular, we recover the behaviour of the depen- 
dence on the cosmology parameters (n m ,n\). However, 
our predictions are somewhat larger than the numerical 
results obtained for low-density universes. This might be 
due to a dependence on Cl m of the parameters S p defined in 
(|Io|). However, as shown by Peacock & Dodds(1996) from 
numerical simulations the relation between the linear and 
non- linear power-spectra for low-density universes is iden- 
tical to the one obtained for the critical density universe 
after one properly takes into account the slower growth of 
the linear density fluctuations. That is the density and the 
size of matter condensations on small non-linear scales, 
which collapsed when Q m ~ 1, is the same as the one 
which would be obtained in a critical density universe as 
seen through £. This suggests that the higher-order cor- 
relation functions and moments should also be the same 
which means that the coefficients S p should not depend 
on O m (the clustering pattern of high-density collapsed 
regions is expected to be roughly the same). 

Moreover, we note that the coefficient S3 which mea- 
sures the skewness of the p.d.f. of the density contrast 5r 
on small scales shows some dispersion with different nu- 
merical simulations. Indeed, while we use S3 ~ 10.7 as 
obtained by Valageas et al.(2000), we note that Colombi 
et al. (1997) find S 3 ~ 10.2 and Munshi et al.(1999) get 
S3 ~ 6. As shown by ( ff9| ) the skewness S Ke> 3 of the conver- 
gence must display the same uncertainty. Thus, it would 
be interesting to measure the coefficients S3 and S Kgt 3 in 
the same simulation. The coefficient S Kf) .3 predicted by 
( ff9| ) increases for larger 9 while numerical results show 
a decline. This decrease is mainly due to the fact that 
for larger angular windows the contribution from linear 
scales grows, which must lead to a smaller skewness 5 Kfl) 3 
since the coefficients S p are smaller in the quasi-linear 
regime (see Colombi et al.1997). The increase with 9 of 
S Ks ,3 (dashed curves) is due to the oscillations of the win- 
dow function W(Vk9) in ([45]). Indeed, for larger angular 
scales 9 the argument T>k9 for k ~ k m , where A 2 (k)/k is 
maximum, reaches values of order unity so that I K is in- 
creasingly sensitive to the oscillations of W(T>k9) (which 
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Fig. 3. The skewness S Kg ,3 of the convergence for a source 
at z s — 1. As in Fig.l we consider three cosmologies: 
open universe, low-density flat model and critical density 
universe, from top downto to bottom. The dashed curve 
corresponds to ( |49| ) and the horizontal dotted line to its 
limit for 9^0 (no smoothing). The slightly lower hori- 
zontal solid line shows the approximation (^). The data 
points are the results of numerical simulations from Jain 
et al.(1999). 

involves the Bessel function J\). As shown by Fig.|| this 
appears for 9 > 1'. Then, the ratio Skg^/S^ given by 
( |64| ) grows as we get further away from the case of a con- 
stant value for the product (wI K ) where the inequality 
( p4| ) would transform to an equality. In this case, the "av- 
erage" provided b y (p9|) is not a very good approximation 
of the expression (|39[) so that the difference between both 
predictions gets larger. However, we can also note that the 
approximation (^9|), which leads to (f30|), is certainly suf- 
ficient in view of the accuracy which can be obtained by 
such methods. The fact that the approximation ( |59| ) shows 
a better agreement with numerical results than the more 
precise expression ( |3^) is a mere coincidence of two effects. 
Indeed, our method may slightly overestimate the skew- 
ness as we do not take into account the fact that the coef- 
ficient S3 gets smaller on larger scales which correspond to 
the mildly non-linear regime. On the other hand, as shown 
in Sect. 3.3 the approximation ( |59| ) leads to a smaller skew- 
ness than the expression ( |39| ) , see (|64|) . This second effect 
makes ( |59| ) to look closer to the numerical points for 9 <; 1' 
but on very small scales the more precise expression ( |39"| ) 
should provide the most accurate predictions. 

Finally, the dependence on the angular scale 9 of the 
skewness S3 shows that for <; 1' the properties of the 
convergence are not perfectly described by the strongly 
non- linear regime (also Hui 1999). Hence one would need 
to describe in details the transition between the highly 
non-linear and linear regimes, so as not to overestimate 
the parameters S p which are larger in the strongly non- 
linear regime (e.g., Colombi et al.1997). Two prescriptions 
have been developed to model this transition: the "ex- 
tended perturbation theory" (Colombi et al.1997) and the 



"hyper-extended perturbation theory" (HEPT) (Scocci- 
marro & Frieman 1999). Both of these empirical models 
use the form of the parameters S p (or of the p— point cor- 
relation functions) in the quasi-linear regime to propose 
an "educated guess" for their behaviour up to the strongly 
non- linear. They have been shown to agree reasonably well 
with numerical simulations. Such a description should be 
extended to all orders to get the p.d.f. P(ng) in the tran- 
sitory regime, which would clearly be quite interesting. 
However, in this article since we focus on small angular 
scales we restrict ourselves to the approximation ( 36 ) with 
constant coefficients S p . Indeed, we shall see below in FigJJ 
and Fig.[| that our predictions for the p.d.f. P{ng) agree 
rather well with the numerical results. This suggests that 
the shape of the p.d.f. P(Kg) is not too sensitive on the 
value of each moment: it rather depends on the behaviour 
of the series of cumulants as a whole. For instance, the pa- 
rameter x s = —l/y s introduced in (A. 2) and (A. 3) which 
governs the high density tail of P{Sr) and the behaviour 
of P{ng) at large Kg, see (pq), is given by: 



0+1 



lim 

p— >oo piip 



(73) 



This small dependence of the p.d.f. of the density con- 
trast 5r (hence of the convergence Kg) on the values of 
the moments (except in the tails of the distribution) was 
already noticed in Valageas et al.(2000). Indeed, in that 
paper the authors found that the p.d.f. P{Sr) obtained 
from various numerical simulations agree fairly well with 
each other (over the range tested by numerical results) 
while the parameters S p display a rather large dispersion, 
as we recalled above. 

4-4- Probability distribution 

Eventually, we can compare our predictions for the p.d.f. 
P(ng) itself with the results of numerical simulations. We 
show in Fig.|| our results at z s = 1 for three cosmologies 
and the angular window 8 = 1'. We can see that our pre- 
dictions (solid lines) agree very well with the numerical 
results (open squares) from Jain et al.(1999). Note that 
at 6 = 1' one might have already expected a small de- 
parture due to the influence of weakly non-linear scales 
where the coefficients S p are slightly different from those 
we use which correspond to the highly non-linear regime. 
In particular, as seen in Fig.||the parameter S K9i 3 we pre- 
dict is not exactly equal to the one which was measured in 
the simulation. As discussed in the previous section, the 
good agreement seen in Fig.[| means that even if our ap- 
proximation for the high-order cumulants (kg) c is not very 
accurate the p.d.f. we obtain can still be a very good ap- 
proximation to the actual P(ng). In particular, it is clear 
from Fig^ that the p.d.f. we obtain is consistent with nu- 
merical simulations for Q m = 0.3, hence the skewness we 
predict, shown in Fig.[| must be consistent with the mea- 
sure from numerical simulations. Thus, the error bars in 
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Fig. 4. The p.d.f. P{ng) of the convergence Kg for a source 
at redshift z s = 1 for three cosmologies and with the angu- 
lar window 9=1'. The solid lines show our prediction ( [42| ) 
using (41) and (p9|). The dashed lines show the gaussian 
which has the same variance (Kg). The dotted lines corre- 
spond to the Edgeworth approximation ([74]). The squares 
show the results of numerical simulations from Jain et 
al.(1999). 

Fig.|| were certainly too small: they represent the disper- 
sion of the estimator used to evaluate S Kg ^ (for instance 
the third moment) in the simulation but other p.d.f. with 
a somewhat different value of S Kgt 3 (e.g., as large as the 
one we obtained) could provide a very good fit the mea- 
sured P{kb). Note also that Fig.[| shows that the generat- 
ing function (p(y) does not strongly depend on the cosmo- 
logical parameters (f2 m , 



Fig. 5. The logarithm of the p.d.f. P(k$) of the conver- 
gence Kg for a source at redshift z s = 1 for three cos- 
mologies and with the angular window 9 = 1'. The solid 
lines show our predictions from the approximation (59) 
as in Fig.||. The dashed lines show the results obtained 
from (|39|). The squares are numerical results from Jain et 
al.(1999). 



We also display in Fig.[| the gaussian (dashed line) 
which would have the same variance (Kg) as our predic- 
tion. This shows that P(kq) cannot be approximated by a 
gaussian. In particular, one can clearly see the asymmetry 
of P(Kg), with a maximum at a small negative value of 
Kg and an extended tail at large positive Kg. Of course, 
this is due to the features of the p. d.f of the density con- 
trast itself. Indeed, on small non-linear scales most of the 
volume is made up of very low density regions ( "voids" ) 
which explains that P(Kg) peaks at a negative value of 
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Kg (corresponding to negative density contrasts along the 
line of sight, see (|27|)). On the other hand, overdense col- 
lapsed matter condensations show an exponential tail at 
large densities which translates into the large-Kg tail of 
P{ne). 

When a p.d.f. shows a small departure from a gaussian 
it is customary to consider the asymptotic Edgeworth ex- 
pansion (e.g., Bernardeau & Kofman 1995). The latter 
is simply obtained by keeping the linear and quadratic 



terms in y within the exponential in (41) while expanding 



the exponential of the higher-order terms (of course for a 
gaussian (pa e (y) — — y 2 /2 and there are no higher-order 
terms). To the first order in we obtain: 



P E {k e ) 



e -«9/(2fft 9 ) 



K 



(74) 



which involves the skewness (through Sa g ^) in the cor- 
rection term. Here H3 is the Hermite polynomial of or- 
der 3 defined by: H 3 (x) = 8x 3 - 12x. The p.d.f. @ is 
shown as the dotted line in Fig.[|. We see that it cap- 
tures the shift of the maximum of P{ng) towards negative 
Kg. However it fails for large Kg where a spurious oscilla- 
tion appears (because of the cubic Hermite polynomial). 
Note that the Edgeworth expansion is only an asymp- 
totic expansion which is useless for <; 0.5. Thus, 
in mildly non-linear regimes one cannot reconstruct the 
p.d.f. from the first few moments. In particular, the re- 
lation ((73]) clearly shows that even to obtain the basic 
features of P(Kg) one needs the complete set of the coef- 
ficients S p . 

We present in Fig.|| a co mpar ison of the two estimates 
of P(Kg) we obtained in Sect. 3.2. First, we note that both 
approximations give results which are very close and show 
a reasonable agreement with the results of numerical sim- 
ulations. As discussed in Sect.|3.3|, we can check that the 
(|5£ 



approximation ( |59| ) leads to a slightly sharper cutoff for 
the exponential tail of P(Kg) than the more precise expres- 
sion ([39"|), see ([33|). However, it is interesting to note that 
the simple approximation (|5^) provides very good results. 
Thus, in view of the inaccuracy introduced by the approxi- 
mation (|36|) and the dispersion of the actual coefficients S p 
(or equivalently of the function h(x) defined in (|ilf)) mea- 
sured from counts-in-cells with different simulations the 



p.d.f. obtained from (59) is probably sufficient. To signifi- 
cantly improve over (59) one would probably need a better 
treatment than 



B_6|) and a better accuracy of the measure 
of the p.d.f. P{5r) or a rigorous theoretical derivation of 
the coefficients S p . The Fig.|5| also clearly shows the asym- 
metry of the p.d.f. P(k$), with a sharp cutoff for negative 
Kg, due to the lower bound K m ; n , and an extended expo- 
nential tail at large positive Kg. 

The fact that the simple approximation (|5^) agrees 
well with the numerical results is quite interesting. Indeed, 
it shows that from the p.d.f. P(Kg) of the convergence Kg 



n 1 1 1 1 1 , , , 1 1 
. iog[ </c| > 1/s /k m , n | ] 




log <K| >VE 



M—\ — 1 — — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — h 



OCDM 




Fig. 6. The dependence on the redshift z s of the source of 
the logarithm of the variances J (k'q) (convergence, dashed 
line) and yj (kg) (normalized convergence, dotted line), 
and of the lower bound |K m i n | (solid line). We consider 
the window 9 — 1' for the three cosmologies. 

we can directly derive the generating function <p(y) which 
describes the non-linear density field at scales defined by 
a comoving wavenumber k ~ 10 Mpc -1 (see Fig.||), that 
is for typical lengths x ~ 0.1 Mpc. Then, from <p(y) we 
can obtain the p.d.f. P(Sr) of the density contrast using 
(pi]). This shows that the p.d.f. P(k$) of the convergence 
is a very efficient probe to obtain the p.d.f. P(5r) of the 
underlying density field. 



5. Dependence on the redshift of the sources 

Finally, we briefly consider the dependence on the redshift 
z s of the source of the weak gravitational lensing effects 
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Fig. 7. The dependence on the redshift z s of the cutoffs 
l^minl (solid line), ng tS ,m (dotted line) and Kg,s (dashed 
line) for the convergence. We consider the window 9 = 1' 
for the three cosmologies. 

we have described in the previous sections. Thus, we show 
in Fig.^| the redshift evolution of the variance yj (kq) and 
its normalized counterpart v/ (M) for the convergence Kg . 
First, we note that the fluctuations yj (ni) increase with 
the redshift z s . Indeed, since the line of sight becomes 
longer as it extends to higher redshifts the distortions 
due to weak gravitational lensing increase. Note however 
that the growth of these fluctuations is rather slow beyond 
z s = 1. This is due to the slow increase of the radial co- 
ordinate x, which appears for instance in (|7l|), because it 
has a finite limit for z — * oo, see (|2p|). Moreover, at large 
redshift the density fluctuations are smaller as structure 
formation is less advanced, which also tends to diminish 



the contributions from high redshifts. On the other hand, 
the "normalized" variance yj (kf) decreases for larger red- 
shift z s . Indeed, the quantity |K m in| (solid lines) increases 
faster with z s than the r.m.s. fluctuation as it is not af- 
fected by the redshift dependence of the power-spectrum, 
see (pf). In particular, as we noticed in Valageas (1999a) 
we obtain from Sect. 3.3: 







and £ ft 



and: 



oo : £ Ke is finite and ^ e oc (1 + z s y 



(75) 



(76) 



Note that a smaller normalized variance yj (kg) implies 
that the p.d.f. is closer to a gaussian near its maximum. 
In particular, at small redshift where \J (k!) JSb |K m in| 
the p.d.f. P(ko) is very different from a gaussian as it 
is strongly affected by the lower bound K m i n . 

Next, we display in Fig.^the evolution with redshift of 
the cutoff which characterizes the exponential tail P(kq). 
First, in agreement with the discussion about Fig|| we 
note that the magnitude |fv m i n | of the lower-bound of the 
convergence increases faster with redshift than the expo- 
nential cutoff which depends on the normalization of the 
power-spectrum at each redshift, as seen for instance in 
( |58|) or (|6ll). We note that at large redshift the cutoff 
Kg yS defined in ( |58|) is much larger than the prediction 
Ke.s,m, defined in (|6l[), of the simple approximation (|59"|). 
This agrees with the discussion of Sect.^, see (|63|), and 
with the behaviour shown in Fig.|| where we can see that 
for very large Kg the p.d.f. P{ng) obtained from ( |59| ) is 
too small as compared to the prediction of the (theoret- 
ically) more accurate relation (|39|). Nevertheless, as seen 
in Fig.^ these two approximations agree very well over a 
large range for Kg. This shows that it is rather difficult 
to measure the actual cutoff of these probability distribu- 
tions (see also Valageas et al.2000 for a discussion for the 
case of the p.d.f. P(<5_r)). 

As we described in Sect .^5| our results also apply to the 
case of a broad redshift distribution of sources n(z s )dz s . 
Moreover, we noticed that since tp(y) is a good approxima- 
tion to (pn e (y) almost all the dependence on the cosmology, 
the source redshift and the angular window, is enclosed in 
the two numbers K m j n and (or £ Ke ) (here we omit the 



bar used in Sectj3.5| to distinguish the case of a "broad" 
redshift distribution of sources). The variance £ K „ of the 
convergence can be directly obtained from the observa- 
tions, hence we are left with only one parameter: K m i n . 
Note that the latter only depends on the cosmological pa- 
rameters (f2 m , Q\, Hq) and on the redshift distribution of 



the sources, see (29): it is independent of the properties 
of the density field (power-spectrum,...). Thus, if these 
cosmological parameters are know (e.g., from CMB data 
and SNela) and the redshift of the sources is well mea- 
sured it can be computed from (p9|). On the other hand, 
if one assumes that the skewness of the density field £3 
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is known to a good accuracy (e.g., from numerical simu- 
lations) then K m ;„ can be obtained from the observations 
using (|6(i|). Bernardeau et al.(1997) advocate the use of 
this method to measure the cosmological parameters, us- 
ing the skewness S3 obtained by perturbative methods 
in the quasi-linear regime relevant to large angular scales 
(9 <; 10'), if the redshift of the sources is well measured. 
However, for the small angular scales studied in this ar- 
ticle (6 55 1') it may be better to obtain K m i n from a fit 
to the full p.d.f. rather than to the skewness only since 
the latter may not be very accu rately constrained (see the 

and Valageas et al.2000). 



discussion in Sect. 4.3, Sect 



4.4 



Nevertheless, the observations of weak gravitational lcns- 
ing effects on small scales should probably be seen as a 
tool to constrain the properties of the density field rather 
than the cosmological parameters. 

6. Conclusion 

In this article we have described a method to obtain an 
estimate of the full p.d.f. of the convergence n. We sum- 
marize below the main new results we have obtained. 

- Extending the method used in Valageas (1999a) we 
show how one can express the moments of the convergence 
in terms of the moments of the density field at all orders, 
working in real-space. This was already studied for the 
moments of order two and three in the non-linear regime 
(e.g., Hui 1999). Since our approximation holds for all or- 
ders it allows us to sum up the series of the cumulants and 
to write the full p.d.f. P(kq) of the convergence in terms of 
the p.d.f. P(5r) of the density contrast on comoving scales 
x ~ 0.1 Mpc (for angular windows 6 55 1'). This provides 
an explicit link between two properties of the density field: 
the counts-in-cells statistics and the convergence of weak 
gravitational lensing distortions. Our results apply to the 
case of a broad redshift distribution of the sources as well 
as to the case where all sources are located at the same 
redshift. 

- Then, we compare our predictions with the results of 
numerical simulations. As pointed out by previous studies 
(e.g., Jain & Seljak 1997) we find that on these angu- 
lar scales (9 <, 1') the non-linear evolution of the den- 
sity field is important. We show that our predictions for 
P(ng) agree very well with the results from N-body simu- 
lations, although there is a non-negligible discrepancy for 
the skewness S Kg ^. In particular, one would need a more 
detailed treatment (i.e. to model the transition between 
the highly non-linear and linear regimes) and more accu- 
rate simulations (the estimate they give for S3 still varies 
from 6 to 10.7) to use the skewness on small angular scales 
to measure the cosmological parameters (the situation is 
easier at large angular scales where one can obtain £13 from 
perturbative methods). On the other hand, it might be 
useful to fit the full p.d.f. rather than the skewness alone. 
Note that on small angular scales one cannot reconstruct 
P(k$) from its first few moments (variance and skewness) 



as the Edgeworth approximation is only an asymptotic 
expansion. Moreover the exponential tail of the p.d.f. is 
significantly different from a gaussian. 

- We have also devised a very simple approximation 
which gives very good results for P(k$) despite its detailed 
properties (location of the singularity at y s , value of the 
skewness) are somewhat different from the prediction of 
the more detailed prescription. It shows that the measure 
of P(ne) provides a direct estimate of P(Sr). Moreover, 
it also means that the p.d.f. P(ng) should obey a specific 
scaling property to a good accuracy. Indeed, to a good 
approximation the p.d.f. P(kg) of the normalized conver- 
gence kg should only depend on one number: the variance 

. Hence the p.d.f. obtained for various angular windows 
9i at a redshift z s \ should nearly superpose onto the p.d.f. 
obtained for another redshift z s i after the rescaling such 
that — £k82 ■ 

- Finally, we have briefly presented some aspects of 
the dependence of gravitational lensing distortions with 
the redshift of the source. 

Using the relations between the properties of the den- 
sity field and the gravitational lensing effects we have de- 
scribed in this article, one could constrain the cosmologi- 
cal scenario of structure formation from observations. For 
instance the ellipticities of images of galaxies can allow 
one to measure the shear 7. Then, one could derive both 
the cosmological parameters (f2 m ,f2A) (e.g., through the 
second and third order moments) and the p.d.f. of the 
density field on non-linear scales. This could provide very 
interesting results, complementary to the usual surveys 
of galaxies which give indirect constraints on the mass 
functions. Since the latter can be related to the statis- 
tics of the counts-in-cells (e.g., Valageas & Schaeffer 1997; 
Valageas et al.2000) this would provide a good check of 
our descriptions of the density field. We shall also present 
in future articles (Valageas 2000; Bernardeau & Valageas 
2000) how one can extend our results to other measures of 
weak gravitational lensing, like the shear or the aperture 
mass. 

However, we must note that in this article we did not 
take into account several effects which might distort the 
p.d.f. P(k$) we obtained. Thus, the non- linear coupling 
between deflecting lenses along the line of sight and the 
higher-order terms beyond the Born approximation lead 
to non-linear terms in the expression of the convergence as 
a function of the density contrast 5 which would slightly 
change the p.d.f. P{ne). These effects have been studied 
in the quasi-linear regime (i.e. 9 10') where they have 
been found to be negligible for low-order moments (e.g., 
Bernardeau et al.1997). A detailed study would be needed 
to estimate their magnitude in the non-linear regime but 
the good agreement of our predictions with the results 
of numerical simulations suggests that these corrections 
should be rather small (and probably below the inaccuracy 
of these simulations). On the other hand, the coupling 
between the sources and the lenses, which we mentionned 
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in Sect. 3.5, is not taken into account in the numerical Appendix A: Density probability distribution 



simulations. This was studied in the quasi-linear regime 
in Bernardeau (1998) who found that the correction is 
negligible for the skewness and the kurtosis if the width 
of the redshift distribution is small enough (Az s ^ 0.15). 
Although we may expect similar results at small angular 
scales it would be interesting to investigate in details this 
effect. 
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Here, we recall the behaviour of P(6r) obtained from sim- 
ple forms for h(x) which are consistent with numerical 
simulations. From very general considerations (Balian & 
Schaeffer 1989) one expects the function (p(y) defined in 
(|l0|) to behave as a power-law for large y: 



y — >• +oo : ip(y)^ay 1 ^ UJ with < ui < 1 , a > 0(A.l) 

and to display a singularity at a small negative value of y: 

y^yf : <p(y) = -a. T(lo s ) (y - y s )-"° (A.2) 

where we neglected less singular terms. From this be- 
haviour of f{y) we have (Balian & Schaeffer 1989) from 



, . , a(l — u>) , , o 
x> 1 : h{x) ~ a s x"^ 1 e~ x/Xs 



(A.3) 



with x s = l/\y s \. Thus, using Jl5| ) we can see that the 
density probability distribution P(5r) shows a power-law 

behaviour from (1 + 8r) ~ £ U up to (1 + Sr) ~ x s £, 

with an exponential cutoff above x s £. At present, there is 
no good theoretical model to obtain the functions h(x) 
and f(y). In particular, Valageas & Schaeffer (1997) and 
Valageas et al.(2000) compared the functions h(x) given 
by numerical simulations with the predictions of a sim- 
ple model (consistent with the usual Press-Schechter pre- 
scription, Press & Schechter 1974) based on the stable- 
clustering ansatz and the spherical collapse dynamics. 
They found that although this model recovers the quali- 
tative change of h(x) with the slope n of the linear power- 
spectrum it strongly disagrees with numerical results on a 
quantitative level. Similar results were reached in Valageas 
(1998) where the full p.d.f. of the density contrast (for 
voids as well as for overdensities) was compared with more 
detailed spherical models and with the adhesion model. 
Thus, in this article we use the scaling function h(x) ob- 
tained from numerical simulations by Valageas et al.(2000) 
for the case of a critical universe with an initial linear 
power-spectrum which is a power-law n = — 2: 



h(x) = 



with: 



a(l-w) x u 



T(lu) (l + bx) c 



eiq>(-x/x s ) 



(A.4) 



a = 1.71 , lj = 0.3 



= 13 , 6 = 5 and c = 0.6 



hence uj s = uj — 1 — c = —1.3 

Indeed, as seen in Figj^ the contributions to weak lens- 
ing distortions come from scales where n ~ —2. However, 
the curvature of the CDM power-spectrum may slightly 
change the parameters S p from the value they would have 
for a pure power-law P{k). Nevertheless, the use of (A.4) 
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has the advantage of the simplicity and should provide a 
reasonable description for all cases of interest, as seen in 



Fig.[4|. The scaling function h(x) shown in (A. 4) defines the 
generating function ip(y) through (|l8|). In particular, one 
obtains (see Gradshteyn & Ryzhik 1965, §9.211, p.1058): 



and 



tp(y) = ab 



l-UJ 



ip [ui — 1, 1 + U) s 



y-y s 



y - y s 



(A.5) 



(A.6) 



-ip[u - l,l + u s ;—r— 



where ip is Kummer's function which can be expressed in 
terms of the difference between two confluent hypergeo- 
metric functions \F\. As emphasized in Valageas (1999a) 
it is better to define tp(y) from h(x) rather than trying to 
use a fit for <p(y) itself. Indeed, from (|Io|) we see that: 



p > 1 



(_i)p-V (p) (o) = s p >o 



(A.7) 



and moreover one can show from ( |r^ ) that the coefficients 
S p must obey: 



S 



p+q dp—q 



(A. 



These constraints are automatically verified if one defines 
f{y) from h(x). If one uses a fit for (p(y) which does not 
obey these constraints one may get negative probabilities 
(since in this case h(x) has to be negative in some range). 
Moreover, the expansions around y = of ip(y) and <p(y) 
must satisfy: 



f{y) = y 



2 



<p{y) 



T_ 

' 2 



(A.9) 



This is important in order to get the right location of the 
maximum of the various p.d.f. we consider as well as the 
right variance. 

Appendix B: Tree-model 

A popular model for the p— point correlation functions in 
the non-linear regime is to consider a tree-model (Schaeffer 
1984; Bernardeau & Schaeffer 1992) where £ p is expressed 
in terms of products of £2 as: 



= E ( 2p Q) Ell^,r J ) 

(«) 



(B.l) 



p-i 



where (a) is a particular tree-topology connecting the 
(p — 1) points without making any loop, Qp°^ is a pa- 
rameter associated with the order of the correlations and 
the topology involved, t a is a particular labelling of the 
topology (a) and the product is made over the (p— 1) links 
between the p points with two-body correlation functions. 



In this case, the coefficients S p defined in ( pL Q| ) are given 
by: 



„ _ \ ' n (q) \ ' Jv V? P Hp-1 &( T i' T 3 



(n) 



(B.2) 



On the other hand, the convergence Kg involves integrals 
along the line-of-sight. In particular, we need the quanti- 
ties Up defined by (see (]3F 



/OC P 
Y[dxi £ P 
-°°i=2 





VC1 



(B.3) 



Then, it is easy to see that the angular functions uj p dis- 
play the same tree-structure (B.l) as the 3D correlation 
functions £ p (see Bernardeau 1995 for a detailed study in 
the quasi- linear regime): 

Wp (Cx, C P ; Qp a) E II <"»(£*> Cj',z) (B.4) 



(a) 



with: 



r°° dk A 2 (k) 

W2(Cl,C 3 5«)=7rjf y ^f 1 Jo (Vk\d ~ C2I) (B.5) 

The cumulants (/tg) c depend on the mean over the disk 
of radius 6 of the quantities cu p , see (|36|), that is on the 
coefficients s p given by: 



(a) 



(a) 



E 



(7r6» a )P 



n P -i w 2 (Ci,c,-) 



(7r6» 2 ) 2 



w 2(Cl,C 2 ) 



(B.6) 



Note the similarity with (B^). However, since the vol- 
ume averages are different (a disk or a sphere) the rela- 
tive weight of the various topologies are slighlty different 
in (B.6) as compared to (B.2), so that s p cannot be ob- 
tained directly from S p . In ([36D we simply made the ap- 
proximation Sp ~ S p which should be quite reasonable. 
For instance, Gaztanaga (1994) obtains S3 ~ 0. 98153 for 
a power-law two-point correlation function £(r) oc r~ 1,8 . 
Note moreover that we can expect the approximation (36) 
to be reasonable for models more general than the tree- 



model (B.l). Furthermore, the inaccuracy introduced by 



the approximation s p ~ S p is much smaller than the un- 
certainty which affects the measure of S p from counts-in- 
cells statistics, as noticed in Sect. 4.3. 



